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Abstract. I review some important results from RHIC experiments. They were obtained in a unique envi- 
ronment for studying QCD bulk matter at temperatures and densities that surpass the limits where hadrons 
can exist as individual entities, raising the quark-gluon degrees of freedom to prominence. These findings 
support the major experimental observations from measuring the bulk properties of particle production, 
particle ratios and chemical freeze-out conditions, elliptic flow followed by hard probes measurements: 
di-jet fragment azimuthal correlations, high-p T hadron suppression, and heavy-flavors probes. I present the 
measurements as a function of collision centrality, energy, system size and for different particle species. 
These results reveal that a dense strongly interacting medium was created in central Au + Au collisions 
at v / snn = 200 GeV: the RHIC discovery. Further, they suggest that this medium is partonic. However, 
the discoveries so far observed at RHIC are far from being understood fully. Accordingly, the focus of the 
experiments has shifted from the discovery phase to the detailed exploration phase of the properties of this 
medium. 

PACS. 25.75.-q - Relativistic heavy-ion collisions, 25.75.Ag - Global features, 25.75.Ld - Collective flow, 
25.75.Gz - Particle correlations, 25.75.Bh - Hard scattering, 25.75.Cj - Heavy quark production, 25.75.Nq 
- Quark-gluon plasma production 



1 Physics motivation and RHIC achievements 

Quantum Chromodynamics (QCD) is considered the fun- 
damental theory for strong interactions. According to it, 
hadronic matter under extreme dense and hot conditions 
must undergo a phase transition [1] to form a Quark Gluon 
Plasma (QGP) in which quarks and gluons no longer are 
confined to the size of a hadron. Recent results from lat- 
tice QCD at finite temperature [3J[3] reveal a rapid in- 
crease in the number of degrees of freedom associated 
with this deconfinement of quarks and gluons from the 
hadronic chains. The transition point is at a temperature 
T « 170 MeV and energy density e sa 1 GcV/fm 3 . Under 
the same conditions, chiral symmetry is restored [2] . There- 
fore, experiments search for signatures of both QGP for- 
mation and the in-medium effects of hadron properties. 
It was proposed that the required high densities could be 
achieved via relativistic heavy ion collisions [3] . 

Under RHIC (Relativistic Heavy Ion Collider) project, 
an accelerator was constructed at Brookhaven National 
Laboratory (BNL) from 1991 to 1999. RHIC was designed 
as a heavy-ion machine, able to support the collision of a 
wide range of nuclei over a large range of energies. Al- 
ready, gold-gold, copper-copper and deuteron-gold colli- 
sions have been attained at energies from ^/snn = 19.6 
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to 200 GeV. Further, a polarized capability was added 
to RHIC allowing transverse and longitudinal polarized 
protons to collide at energies from ^/snn = 200 to 410 
GeV. By now, RHIC has made a major physics discovery, 
namely the creation in high-energy central gold-gold col- 
lisions of a new form of matter; dense and strongly inter- 
acting, called the strongly coupled quark-gluon plasma, 
or sQGP. This finding was rated the top physics story 
of 2005 and the four experiments at RHIC, BRAHMS, 
PHENIX, PHOBOS, and STAR, published, in white pa- 
pers, evidence of the existence of this new form of matter 
[6]. RHIC accelerator also met and surpassed its specifi- 
cations; namely, it attained its energy goals, and exceeded 
the heavy- ion luminosity goals by factor of 2, and polar- 
ized proton luminosity by a factor of 5. Table [T] summa- 
rized the achievements of RHIC over the last eight years. 

In this article, I highlight the recent RHIC results un- 
derlying the major experimental observations, i.e., "dis- 
coveries". I start with discussions of bulk particle pro- 
duction and the initial conditions, followed by particle 
ratios and chemical freeze-out conditions. I then present 
our measurements of elliptic flow that are an indirect sig- 
nature of the existence of partonic matter, followed by 
hard-probe measurements: di-jet fragment azimuthal cor- 
relations, high-p T hadron suppression, and heavy-flavors 
probes. These measurements of hard probes afford direct 
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Fig. 1. Panel a) The measured distribution of charged-hadron as a function of transverse momentum p T for the 15% most 
central Au + Au collisions at -^/snn = 200 GeV. Panel b) The measured scaled pseudorapidity density dNch/d7//(^N P art) 
for | r] [< 1 in central Au+Au collisions at the AGS, in central Au + Au and Cu + Cu collisions at the RHIC and in cen- 
tral Pb + Pb collisions at the SPS [5H5l llO|[TT] . The star symbols denote the extrapolation for dN c h/dr7/(iN par t} in central 
Pb + Pb collisions at 5.5 TeV (LHC). The continuous curve corresponds to the logarithmic quadratic fit of the RHIC data: 
f^ A = 3.09 - 1.061n(^/S^7) + 0.22(ln(^s7^7)) 2 . The dashed line represents the logarithmic linear fit of the AGS, SPS and RHIC 
data: f£ A = 0.4749 + 0.77 ln( v /s^7). The error bars show the systematic errors. The Color Glass Condensate model (CGC 
model) prediction for LHC is also illustrated as a red box [P2] . 



signatures of highly interacting dense matter created in 
central Au + Au collisions at v / Snn = 200 GeV. The mea- 
surements, for different particle species, are given as a 
function of collision centrality, energy, and system sizes. 

2 Global properties of hadron production 

Measurements of charged-particle multiplicity and trans- 
verse energy distributions, and of azimuthal anisotropies 
in heavy-ion collisions afford information on the initial en- 
ergy density and the entropy production during the sys- 
tem's evolution, and are sensitive to a variety of physics 
processes responsible for generating multiparticles. This 
knowledge is important for constraining model predictions 
and indispensable for understanding and estimating the 
accuracy of the more detailed measurements of, for ex- 
ample, jet or quarkonia production. Figure [Ij, shows the 
measured distribution of charged-hadrons produced in the 
15% most central Au + Au collisions at ^/snn = 200 GeV [7] 
The distribution illustrates a "bulk" and "tail" that are 
respectively related to "soft" and "hard" parton-parton 
scattering. We clearly observed that the "bulk" consti- 
tutes the dominant part of charged-hadron production. 
However, there is no clear separation between "soft" and 
"hard" processes. For this reason, an analysis undertaken 
as a function of centrality, transverse momentum, energy, 
and system sizes is required for better understanding par- 
ticle production. 



2.1 Charged-particle density distributions and initial 
conditions 

One important observable in heavy-ion interactions is the 
number of charged-particles produced in a collision, called 
the charged-particle pseudorapidity density dNch/dry. This 
number (dN c h/df?) is proportional to the entropy density 
at freeze-out and, since entropy cannot be destroyed (even 
in non-equilibrium systems), pseudorapidity density pro- 
vides information on the initial state of parton density, 
and any further entropy produced during subsequent evo- 
lution [5]. 

Fig. [1}d shows the charged particle densities near the 
mid-rapidity region, dN ch /d?7|^| <1 /(Np art /2), where (N part ) 
is the average number of nucleon participant pairs, for 
Au + Au (RHIC: 19.6, 62.4, 130 and 200 GeV), Cu + Cu 
(RHIC: 22.4, 62.4 and 200 GeV), Pb + Pb (SPS: 8.83, 
12.28 and 17.26 GeV) and Au + Au (AGS: 2.63, 3.28, 
3.84, 4.29, 4.86 GeV); they present the 6% most cen- 
tral collisions as a function of center-of-mass collision cn- 
er gy (v^nn) [51I WTU1HT] . The results in this figure (see 
Fig. Qi>) suggest that the particle density rises approxi- 
mately logarithmically with energy. Comparing the find- 
ing for Cu + Cu and Au + Au indicates that, for the most 
central events in symmetric nucleus-nucleus collisions, the 
particle density per nucleon participant-pair does not de- 
pend on the size of the two colliding nuclei but only on 
the collision's energy p~0j[TT]. 
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Fig. 2. The measured distributions of pseudorapidity density in Au + Au, Cu + Cu, d + Au, and p + p collisions at RHIC 
energies 8,9, 10; TT|. The dN c h/d?7 distributions for Au + Au, Cu + Cu and d + Au collisions are plotted as a function of collision 
centrality. Typical systematic errors (90% C.L.) are represented as bands for selected centrality bins. Statistical errors are 
negligible. 



Based on the data presented in the Fig. [TJd, we can 
establish two types of fits to predict (extrapolate) the re- 
sults of dN c h/d?y at the mid-rapidity region in Pb + Pb 
collisions at the LHC energy 5.5 TeV: 

1. The logarithmic quadratic fit of the RHIC data is shown 
as the solid curve in Fig. [TJd and its fit function is: 



f Q 

AA 



3.11 - l.07ln( 1l 



0.23(ln( 



(1) 



this formulation allows us to extrapolate the scaled 
density per nucleon participant pair for central Pb + Pb 
collisions at LHC energy (5.5 TeV): 



1 



x — — (Pb 



(N part /2) dr? 



Pb at 5.5 TeV) = 10.95 



Using a Glauber model calculation for the 6% most 
central Pb + Pb collisions at 5.5 TeV (the total inelas- 
tic cross section used in the Glauber model calculation 
is a NN — 64 mb), we obtained the value of (N part )= 
365.5, from which we can deduce the value of the un- 
sealed pseudorapidity density at the mid-rapidity re- 
gion: 



dN ch 
dry 



(central : - 6%, Pb + Pb at 5.5 TeV) = 2002 



2. The logarithmic linear fit of the data from the AGS, 
SPS and RHIC is shown by the dashed line in Fig. [TJj, 
and has the fit function: 



f AA = -0-33 + 0.75 ln( A 



(2) 



With this fit equation, we extrapolated the scaled par- 
ticle density per nucleon participant pair at the mid- 
rapidity region to LHC energy, ^/snn = 5.5 TeV: 



1 



x — — (Pb 



(Npart/2) dry 



Pb at 5.5 TeV) = 6.1 



For (Np ar t) 
calculation 
collisions at 5.5 TeV, this yielded 



= 365.5 obtained via the Glauber model 
using <r NN = 64 mb) for 6% central Pb + Pb 



dN. 



ch 



drj 



(central : - 6%, Pb + Pb at 5.5 TeV) = 1120 



The prediction of the Color Glass Condensate model (CGC 
model) pj] for dN ch /cIt? 1 1 | < i / (N part /2) in Pb + Pb colli- 
sions at LHC energy, -^/snn = 5.5 TeV, agrees very well 
with the results of extrapolation obtained with the loga- 
rithmic quadratic fit, as shown in Fig. [TJa. 

Fig. [H plots the measured dN c h/dry of primary charged 
particles over a broad range of pseudorapidity, \rj\ < 5.4, 
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Fig. 3. Panel a) compares the experimental data from different particle multiplicity ratios obtained in Au + Au collisions 
at RHIC energies at ^/snn = 130 and 200 GeV with calculations from the thermal model [19] • Panel b) shows the ratios of 
antiparticles to particles as a function of rapidity |20| . The errors bars depict the statistical errors and the caps indicate the 
combined statistical errors. Panels c) and d), respectively, illustrate the thermal model fit baryon chemical potential (fJ, B ) and 
the temperature chemical freeze-out (T c h) as function of collisions centrality ((N part }) [21]. fj, B and T c h are obtained from 
Au + Au, Cu + Cu and p + p collisions at RHIC energies and for different particle species. 



for Au + Au [Eli] and Cu + Cu flTHITT] collisions under 
a variety of collision ccntralities and RHIC energies. The 
d + Au p~3t[14:. 15J and p + p [13] data at RHIC energies 
also are shown. The Cu + Cu, Au + Au d + Au and p + p 
data at all energies were obtained with the same detector 
setup in the PHOBOS experiment [J6]. This is optimal be- 
cause common systematic errors cancel each other out in 
the ratio. With this configuration, we were able to exam- 
ine comprehensively particle production in Cu + Cu and 
Au + Au collisions for the same number of nucleon par- 
ticipant pairs, for the same fraction of total inelastic cross 
sections, and for the same geometry in both systems |llj . 
Both the height and width of the dN c h/dr; distributions 
increase as a function of energy in both systems. Com- 
paring the results revealed an interesting feature, namely 
that the best agreement of dN c h/(N part /2) distributions 
over the full coverage, \rj\ < 5.4, for central and for pe- 
ripheral Cu + Cu and Au + Au collisions is obtained for 
centrality bins selected to yield a similar geometry, i.e., 
a comparable value of N part /2A (where A is the atomic 
number) in both systems [TTj . 

One question to be asked is whether mechanisms in 
the limiting fragmentation region (i.e. high \r)\) are dis- 
tinct from those at mid-rapidity (jr)\ < 1). From the results 



shown in Fig. [5] for Au + Au, Cu + Cu, d + Au, and p + p 
collisions, there is no obvious evidence for two separate re- 
gions at any of the RHIC energies [TO]. However, to date 
no such anomalies have been noted, at the AGS, SPS, or 
RHIC. So far, all results on particle pseudorapidity densi- 
ties point to a rather smooth evolution in centrality and 
y's NN . This of course, does not necessarily imply the ab- 
sence of a phase transition, but rather, might indicate the 
insensitivity of these observables to the early phase of the 
collision and/or might suggest a second order phase tran- 
sition (or a cross-over) jT7] . We note that charged-particle 
production mechanisms for dN c h/d?7 shown in Fig. [5] are 
mostly obtained from the "soft" processes (see Fig. [JJi). 

2.2 Particle ratios and chemical freeze-out conditions 

The chemical freeze-out point is the stage in the evolu- 
tion of the hadronic system when inelastic collisions cease 
and the relative particle ratios become fixed; this point, 
is defined by the temperature chemical freeze-out, T ch , 
and also by the baryon chemical potential, /i B . These pa- 
rameters, T ch , and /i B determine the particle composition 
of the hadronic final state. After chemical freeze-out, the 
particle composition inside the fireball is fixed, but elastic 
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collisions keep the system intact until the final, thermal 
freeze-out. At this stage the momentum distributions of 
particles are final and no longer change. Therefore, the 
transverse momentum spectra determine the parameters 
of the thermal freeze-out. Thereafter, statistical interpre- 
tation of particle production becomes an appropriate ap- 
proach for evaluating heavy-ion collisions at high energies; 
because large multiplicities of hadrons are created. We can 
assume that the nuclear matter created in these collisions 
form an ideal gas that can be characterized by a grand- 
canonical ensemble. Using thermodynamic concepts to de- 
scribe multiparticles production has a long history [TB] . 
The concept of a temperature applies only to systems in 
at least local thermal equilibrium. The assumption of a 
locally thermalized source in chemical equilibrium can be 
tested by applying statistical thermal models to describe 
the ratios of various emitted particles. 

Figure [3^ compares the RHIC's experimental parti- 
cle ratios and statistical thermal model calculations for 
Au + Au collisions at 130 and 200 GeV [15]; note that 
the measurements were taken in the mid-rapidity region 



| r] | < 1. They demonstrate quantitatively the high degree 
of equilibration achieved for hadron production in central 
Au+Au collisions at RHIC energies. Values were obtained 
for (T ch , n B ) of (174 ±7, 46±5) and (177±7, 29±6) at 
130 and 200 GeV, respectively, The model matches well 
with these results, and there is no indication of a signifi- 
cant change in T ch at the two energies, 130 and 200 GeV. 
There is a drop in fi B from ~ 46 MeV, at ^/snn = 130 
GeV to 29 MeV at y/sm = 200 GeV. We note that ratios 
involving multi-strange baryons are well reproduced, as is 
the ratio. Even relatively wide resonances such as 

the K*'s fit well with the picture of a chemical freeze-out. 

The BRAHMS collaboration reported interesting re- 
sults on particle ratios as a function of rapidity [20j . Fig- 
ure [3p shows the tt~/tt + , K~/K + and p/p ratios as a 
function of this parameter. We observe that the ir~ 
ratio is consistent with unity over the considered rapidity 
range, while the K~ / ' K + ratios drops almost by 30% at y 
= 3 from its mid-rapidity value, and the p/p ratio by al- 
most 70%. Fig. and c show thermal model fit tempera- 
ture and baryon potential at chemical freeze-out as a func- 
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tion of collision centrality ((N part )) in Au + Au, Cu + Cu 
and p + p collisions at ^/snn = 200 GeV [21] . Tempera- 
ture seems constant as a function of (N part ) and appar- 
ently independent of system size (Au + Au and Cu + Cu 
systems). The baryon chemical potential exhibits small 
variations as a function of N part , consistent with of be- 
ing independent of system size in Au + Au and Cu + Cu 
collisions. 



3 Evidence of partonic matter at RHIC 
3.1 Anisotropic flow of hadrons 

The anisotropic flow of hadrons has been studied exten- 
sively in nucleus- nucleus collisions at the SPS and RHIC 
as a function of pseudorapidity, centrality, transverse mo- 
mentum and energy [6T f2"!?l[2"9"] . In non-central collisions of 
heavy ions at high energy, the configuration space anisotropy 
is converted into a momentum space anisotropy. The dy- 
namics of the collision determine the degree of this trans- 
formation. For a symmetric system, like Au + Au, the 
second Fourier expansion is a good parameterization of 
the anisotropy. At RHIC, a strong anisotropic flow (v 2 ) 
was observed for all hadrons measured suggesting that a 
strongly interacting system was created in the collisions. 

Hydrodynamic models that assume the formation of 
a QGP were used to model the behavior of the medium 
so created [23]. Comparison of the data obtained at mid- 
rapidity region presented by the STAR, PHOBOS, and 
PHENIX collaborations [6)129] to the hydrodynamic model 
[23 affords strong evidence that the originated medium 
behaves as a near ideal fluid. 

Fig. H] shows the anisotropic flow distributions, V2, for 
identified hadrons it, Kg, p, A, H, 51 and cf> [24 . Calcula- 
tions from a hydrodynamic model are depicted in dashed 



lines, while the dot-dashed lines represent the results of fit 
for quark number scaling n q (for mesons, n q = 2; and, for 
baryons: n q = 3). In the lower p T region, p T < 2 GeV/c, 
the value of V2 is inversely related to the mass of the 
hadron, that is, characteristic of hydrodynamic collective 
motion in operation, see Fig [4ji and b. At the interme- 
diate p T region, the dependence is different. Instead of a 
mass dependence, there seems to be a hadron type de- 
pendence, Fig. Hfb). An interesting result concerns the (j> 
meson, whose mass is close to that of p and A. The V2 
distribution illustrated as circles in Fig.[4}j, is the same as 
the V2 for n and CI. It is known that the 4> meson does not 
participate as strongly as others do in hadronic interac- 
tions, nor can </> mesons be formed via the coalescence- like 
K + + K~ process in high energy collisions [55]. Hence, 
the strong V2 we recorded must have developed before 
hadronization. This observation, together with the f2 V2 
results offer clear evidence for partonic collectivity |26) . 
To demonstrate the scaling properties of v 2 , the follow- 
ing transformation have been performed. The measured 
V2 was scaled by the number of valence quarks in a given 
hadron. For mesons and baryons, respectively, they are 
n q = 2 and 3. The p T also was scaled with the same n q . 
The results are shown in Fig. [4J;. To include the effect of 
the collective motion, the p T was transfered to the trans- 
verse kinetic energy KEt = mx— mass = •v/p^ — m 2 — m 
scaled by n q , where m mass of the particle. The outcome 
is depicted in Fig. [JJl. All of the hadron v 2 scaled nicely 
up to (mx — mass)/n q ~ 1.2. Hydrodynamic calculations 
[27] , shown as dashed- lines in the plot, also are scaled in 
the low p T region. At high-p T , the values of V2 appear to 
fall. These observations about scaling reveal that before 
hadronization, quarks already have acquired the collective 
motion V2, and that when they coalesce, the V2 is passed 
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Au + Au and d + Au collisions. Panel b) compares Raa with RdA for n° and r\ from central Au + Au and d + Au collisions. 
Panel c) compares the Raa of Direct 7 with Raa of 7t°, r\ and in central Au + Au collisions. The error bars correspond to 
the statistical errors. For clarity, the systematic errors are shown as vertical bands. 



to newly formed hadrons. Again this observation verifies 
partonic collectivity. 



4 Signatures of dense matter at RHIC 

4.1 Di-jet fragment azimuthal correlations 

The study of high transverse momentum hadron produc- 
tion in heavy ion collisions at RHIC provides an experi- 
mental probe of the QCD matter in the most dense stage 
of the collisions, wherein quark-gluon deconfincmcnt is 



likely to occur [28] . In particular, two- hadron azimuthal 
correlations support the assessment of back-to-back, hard- 
scattered partons that propagate in the medium before 
fragmenting into jets of hadrons, thereby serving as a to- 
mography probe of the medium. 

Fig. [5] shows measurements of two-particle correlations 
[30] given a trigger particle plunder the condition 4< 
p^<6 GeV/c, along with an associated particle with p^: 
1) for 2 < p(, a) < p^GeV/c as shown in Fig. EK, and 2) 
for 0.15 < p^< 2 GeV/c as shown in Fig. [5b. The two- 
particle correlation are presented as a function of the dif- 
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ference between the azimuthal angles of the two particles 
(A</>) produced in Au + Au, d + Au and p + p collisions 
at v / Snn = 200 GcV. These measurements were made at 
the mid-rapidity region < |Ar/| < 1.4 

Fig. [5Jl demonstrates that the trigger-side correlation 
peak in central Au + Au collisions apparently is the same 
as that measured in p + p and d + Au collisions but the 
away-side jet correlation in Au + Au has vanished. This 
observation is consistent with a large energy loss in the 
medium causing it to be opaque to the propagation of high 
momentum partons. However, it is noticeable in Fig. [3}}, 
which has a wide pM, 0.15 < p^< 2 GeV/c, that the 
away side jets have not really disappear, but simply lost 
energy so that the away-side correlation peak has became 
much wider than that of the p + p and d + Au collisions. 
These collective measurements, from p + p, d + Au and 
Au + Au of two-particle correlations, point to the creation 
of a dense medium in central Au + Au collisions at ^/snn 
= 200 GeV. 



4.2 High-p T hadron suppression: jet-suppression 

In heavy ion collisions (from AGS to RHIC energies) hadron 
production at the mid-rapidity region (|y| < 1.5) rises with 
increasing collision energy. At RHIC, the central zone is 
almost baryon free |31j . Particle production is large and 
dominated by pair production, and the energy density 
seems to exceed significantly that required for QGP for- 
mation [BJ. 

Recently, RHIC experiments revealed suppression of 
the high transverse momentum component, p T of hadron 
spectra at the mid-rapidity region in central Au+Au col- 
lisions compared to scaled momentum spectra from p+p 
collisions at the same energy, v / snn = 200 GeV [BJ . This ef- 
fect, originally proposed by Bjorken, Gyulassy and others 
[32] rests on the expectation of a large energy loss of high 
momentum partons scattered in the initial stages of colli- 
sions in a medium with a high density of free color charges. 
According to QCD theory, colored objects may lose en- 
ergy by the bremsstrahlung radiation of gluons [33J . Such 
a mechanism would strongly degrade the energy of leading 
partons, reflected in the reduced transverse momentum of 
leading particles in the jets emerging after fragmentation 
into hadrons. The STAR experiment established that the 
topology of high-p T hadron emission is consistent with jet 
emission, so that jet-suppression is a valid concept. 

Fig.[BJshows compilation of data using Refs. [34] for the 
nuclear modification factors measured for different colli- 
sions systems at RHIC energies. The nuclear modification 
factor is defined as: 



Raa(pt) 



yield per A + A collisions 
Nbin x (yield per p + p collisions) 

d 2 N A+A /d PT d?7 



N bin d 2 N p+p /dpTd?7 

It involves scaling measured distributions of nucleon-nucleon 
transverse momentum by the number of expected incoher- 
ent binary collisions, N^n [35] . In the absence of any mod- 
ifications due to the 'embedding' of elementary collisions 



in a nuclear collision, we expect Raa = 1 at high-p T . At 
low p T , where particle production follows a scaling with 
the number of participants, the above definition of Raa 
leads to Raa < 1 for p T < 2 GeV/c. 

In d + Au collisions at RHIC energy ( v / Snn = 200 
GeV) as shown in Fig. [BJi, the RdA for charged hadrons 
(h + + h~)/2 is enhanced. This enhancement was expected 
as a consequence of Cronin enhancements, now under- 
stood as an initial state effect [36] which also is seen in pA 
collisions. The Cronin effect is associated with the initial 
multiple scattering of high momentum partons. In con- 
trast, RHIC experiments discovered a factor of 4-5 sup- 
pression, in central Au + Au (at ^/snn = 200 GeV) . In 
the same context, the RdAu of tt° and r) in central d + Au 
collisions exhibits no suppression of high-p T , in contrast 
with the Raa of ir° and T) in central Au + Au collisions 
that is suppressed, as shown in Fig. [Bb. At p T ^ 4 GeV/c, 
we find a ratio RdAu/R-AA ~ 5. Indeed, the RdA distri- 
bution shows the Cronin type enhancement observed at 
lower energies as in Pb + Pb collisions at 17.3 GeV/c [37J. 

Fig. [Bfc summarizes the present status of Raa for direct 
photons, 7T° and rj in central Au + Au collisions at ^/snn = 
200 GeV. The Raa for direct photons are not suppressed 
because they do not interact strongly with the medium. 
The Raa for both tt° and rj's mesons exhibit the same 
suppression relative to the point-like scaled p + p data by 
a factor of ~ 5 which appears to be constant for p T > 
4 GeV/c while the rj mass is much larger than that of 
7T°. This observation combined with Fig. [4] is clearly an 
indication of partonic nature of suppression. 

This data of Raa were described by theoretical calcula- 
tions of parton energy loss in the matter created in Au+Au 
collisions [38]. From these theoretical frameworks, we have 
learned that the gluon density dN g /dy must be approx- 
imately 1000 39J, and the energy density of the matter 
created in the most central collisions must be approxi- 
mately 15 GeV/fm 3 to account for the large suppression 
observed in the data 00] . 



5 Heavy-flavors as probes for dense medium 

As discussed above, RHIC has produced in high-energy 
central Au + Au collisions at ^/snn = 200 GeV a new form 
of matter that is dense, strongly interacting and can be 
partonic, called sQGP. This new state of strongly inter- 
acting matter can closely resemble a "perfect" fluid, with 
very low viscosity and high opacity explained by the equa- 
tions of hydrodynamics [53J . The research focus now has 
shifted from this initial discovery to a detailed exploration 
of partonic matter. Particles carrying a heavy flavor, i.e., 
charm or beauty quarks, are powerful tool for studying 
the properties of the hot, dense medium created in high- 
energy nuclear collisions: they are generated early in the 
reaction, and subsequently diffuse in the putative sQGP. 

The study of heavy flavor in relativistic nuclear colli- 
sions follows two different approaches: 1) The direct re- 
construction of the heavy flavor meson, and 2) the iden- 
tification of electrons and muons from semi-leptonic de- 
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Fig. 7. Compilation of data on Au + Au collisions at v / snn = 200 GeV from Refs. [34J[41]. Panel a) represents the nuclear 
modification factor of heavy-flavor electrons Raa compared with the Raa of 7r° and n ± in central Au + Au collisions. Panel b) 
considers the anisotropic flow of heavy-flavor electrons v^ F with that of V2 of ty° and ir^ in minimum-bias Au + Au collisions 
-^/snn = 200 GeV. The error bars correspond to the statistical errors. For clarity, the systematic errors are shown as vertical 
boxes and bands. 



cays of such mesons. The PHENIX and STAR experi- 
ments at RHIC explore one or both of these methods. The 
STAR collaboration is directly reconstructing heavy fla- 
vor mesons the decay channel D° — > K~ + tt + in d + Au 
and Au + Au collisions. Employing semi-leptonic decays of 
heavy flavor mesons, such as D° — > e + + K + + v over a 
broad p T range, more efficiently measures charm and bot- 
tom quark production and overcomes the limitation im- 
posed by directly reconstructing such mesons. 

At RHIC, two methods are utilized to measure heavy 
flavor production via semi-leptonic decays: 1) Identifica- 
tion of electrons from the decays of D and B— mesons, 
and 2) identification of muons from D— meson decays. 
The PHENIX and STAR experiments pursue the anal- 
ysis of non-photonic electrons (from heavy flavor) [4T1I42] . 
Electron identification in PHENIX largely is based on us- 
ing the Ring Imaging Chernkov detector (RICH) in con- 
junction with a highly granular EMC (electromagnetic 
calorimeter). Their momentum is derived from the cur- 
vature of the track (due to a magnetic field up to 1.15 
Tesla) reconstructed from drift and pad chambers. STAR 



identifies of electrons is done using information on dE/dx 
and momentum gleaned from the Time Projection Cham- 
ber (TPC) and the Time of Flight (ToF) data for the 
low to moderate p T (p T < 4-5 GeV/c) electrons. The bar- 
rel electromagnetic calorimeter's (EMCAL's) data is used 
for moderate to high-p T (p T > 1.5 GeV/c) electrons [42]. A 
major difficulty in both collaborations' electron analyzes 
is the fact that there are many sources of electrons other 
than the semi-leptonic decays of heavy flavor mesons. The 
main sources of background are photon conversion in the 
detector material (less significant in PHENIX that has 
a lesser amount of material than STAR), and ir° and rj 
Dalitz decays. Other sources of background such as u, <fi 
and p decays also must be taken into account. The back- 
ground sources usually are called photonic electrons. 

Fig.[7]shows compilation of data using Refs. :34l l4Tll42~j 
of the nuclear modification factors and elliptic flow mea- 
sured for Au + Au collisions systems at ^/snn =200 GeV. 
Fig. [7^ compares the nuclear modification factor Raa of 
heavy flavor electrons to n° data, and to tt^ data obtained 
from central Au + Au collisions at v / snn = 200 GeV. We 
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observe clear suppression of heavy flavor electrons in cen- 
tral events in high-p T . For p T > 4 GeV/c, the of 
heavy flavors is surprisingly similar to that for ir° and 
tt^~ . Fig. [7b illustrates the distribution of the anisotropy 
of heavy flavors electrons (v^ F ) as a function the parti- 
cles' transverse momentum, p T , in minimum bias Au + Au 
collisions at ^/snn = 200 GeV. Collective behavior is ap- 
parent in heavy- flavor electrons (v^ F > 0); but however, 
it is lower than v 2 of ir° for p T > 2 GeV/c. 

This data presented in Fig [7] indicates that heavy fla- 
vors strongly coupled with the medium. The observed sup- 
pression at high-p T > 4 GeV/c can be explained in terms 
energy loss from the heavy-flavors into the medium. Ref. 
[4Tj offers different theoretical predictions for heavy-flavor 
electrons in central Au + Au collisions, considering dif- 
ferent energy loss mechanisms. However, I note that in 
Fig. [7b, the data extend only to around p T « 4 GeV/c, 
and there are large statistical errors in the v^ F . This im- 
plies that we can reach a comprehensive conclusion about 
R^4 and v^ only when we have more data on Au + Au 
collisions at ^/snn = 200 GeV is needed and an upgrade of 
PHENIX and STAR detectors is required. Furthermore, 
we need measurements of the particles identified as de- 
cay products from charm or beauty flavors using the dis- 
placement of their trajectories from the collision vertex 
(measuring the distance of closest approach, DC A). For- 
tunately, both detectors PHENIX and STAR already were 
upgrade by building silicon vertex trackers for heavy fla- 
vors gui. 



6 Conclusions 

I discussed measurements from RHIC collisions on bulk 
particle production, the initial conditions of the collisions, 
particle ratios and chemical freeze-out conditions, flow 
anisotropy of particles, di-jet fragment azimuthal correla- 
tions, high-pi hadron suppression and heavy-flavors probes 
as a function of collision centrality, energy, system sizes 
and for different particle species. The measurements sug- 
gest that RHIC discovered a new state of matter in central 
Au + Au collisions at ^/snn = 200 GeV. This new form of 
matter is a hot, dense, strongly interacting and is called 
"sQGP" . Whilst this material, and its discovery, is far be- 
ing fully understood, the focus of the RHIC experiments 
has shifted from the discovery phase to a detailed explo- 
ration phase of the properties of the medium. The hope is 
that Pb + Pb collisions at LHC energies (LHC energy is 
approximately 30 times higher than RHIC's maximum en- 
ergy: LHC-energy/RHIC-energy = 5500/200 = 27.5) will 
create a medium with an initial temperature high enough 
so that the interaction weakens sufficiently to generate a 
a gaseous QGP. 
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Table 1. RHIC operating modes and total integrated lumi- 
nosity delivered to all experiments for each run [5]. 
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